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ABSTRACT 



We critically re-examine our earlier derivation of the effective low energy action for QCD in 
4 dimensions with chiral fields transforming non-trivially under both color and flavor, using the 
method of anomaly integration. We find several changes with respect to our previous results, 
leading to much more compact expressions, and making it easier to compare with results of other 
approaches to the same problem. With the amended effective action, we find that there are 
no stable soliton solutions. In the context of the quark soliton program, we interpret this as 
an indication that the full low-energy effective action must include additional terms, reflecting 
possible modifications at short distances and/or the non-trivial structure of the gauge fields in 
the vacuum, such as (-F^„) ^ 0. Such terms are absent in the formalism based on anomaly 
integration. 



★ Research supported by the Israel Science Foundation administered by the Israel Academy 

of Sciences and Humanities. 
| e-mail: fnfrishm@wicc.weizmann.ac.il 
\ e-mail: ftami@wicc.weizmann.ac.il 
§ e-mail: marek@vm.tau.ac.il 



1 



1. Introduction and Motivation 



In the recent years some progress have been made towards establishing the 
connection between the phenomenologically successful non-relativistic constituent- 
quark model (NRQM), and QCD's fundamental degrees of freedom. Kaplan' 11 
proposed a physical picture combining some of the features of the chiral quark 
model' 21 and the skyrmion' 3 51 . It was postulated that at distances smaller than 
the confinement scale but large enough to allow for nonperturbative phenomena 
the effective dynamics of QCD is described by chiral dynamics of a bosonic field 
which takes values in U(N c xNf). This effective theory admits classical soliton 
solutions. Assuming that they are stable and may be quantized semiclassically, 
one then finds that these solitons are extended objects with spin 1/2, and that 
they belong to the fundamental representation of color and flavor. Their mass is of 
order Aqcd and radius of order 1/ Aqct>' 61 It is very tempting to identify them as 
the constituent quarks. Thus the constituent quarks in this model are "skyrmions" 
in color space. 

It turns out that in 2 space-time dimensions this picture is exact . Thanks 
to exact non-abelian bosonization one can rewrite the action of QCD2 in terms of 
purely bosonic variables, which are chiral fields G U(N c xNf). It is then straight- 
forward to demonstrate that the only non-trivial static solutions of the classical 
equations of motion are those which contain either a soliton and an anti-soliton or 
iV c solitons. The solitons transform under both flavor and color, yet their bound 
states are color singlets and have the quantum numbers of baryons and mesons. 

In four space-time dimensions there is no exact bosonization, and therefore 
any attempt at derivation of a similar picture in four dimensions must rely on 
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certain approximations. In our previous work we have derived the approximate 
low-energy effective chiral lagrangian with target space in U(N c xNf) using the 

[9—12] 

approach based on integration of the anomaly equations . Equivalent results 
were independently obtained in ref. 13, using an a priori different approach. 

The purpose of the current work is to critically re-examine the results of ref. 8, 
with particular emphasis on the question whether the action we have derived can 
support stable, time independent classical solutions. The existence of such solu- 
tions appears to us to be a necessary condition for establishing the physical picture 
in which constituent quarks are solitons of a low-energy effective action of QCD in 
four dimensions. 

We find several changes with respect to our previous results, leading to much 
more compact expressions, and making it easier to compare with results of other 
approaches to the same problem, and in particular with the action proposed by 
Kaplan 111 . 

With the amended effective action, we find that there are no stable soliton 
solutions. In the context of the quark soliton program, we interpret this as an 
indication that the full low-energy effective action must include additional terms, 
reflecting possible modifications at short distances and/or the non-trivial structure 
of the gauge fields in the vacuum, such as (i 7 ^) ^ 0. Such terms are absent in the 
formalism based on anomaly integration. 

The layout of the paper is as follows. In Section 2 we rederive the low energy 
classical action, resulting in much more compact final expressions. In Section 3 we 
examine the necessary conditions for existence of stable, time independent classical 
solutions with finite energy, and conclude that in order for such solutions to exist, 
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the effective action must contain additional terms which are not present in our 
derivation. Section 3 is devoted to discussion and interpretation of the results. 

2. Derivation of the effective action 
We follow the conventions and notation of ref. 8. 

The variation of the determinant under the axial transformation of the external 
fields, including terms up to zeroth power of the cutoff A, is given by 



where = - d v V^ + % [V^ V v \ + % [A^ A v \ = \{L^ + R^), D^A V = d^A v + 
% [V^A V ] and = \(R^ V - L^) = (D^ - D V A^). (cf. eq. (4.1) of ref. 8). 



is the external source coupled to the vector current of the fermions, while A 
is the external source coupled to the axial current of the fermions. 

An action which has equation (1) as its variation is given by 



—i 




(1) 





Tr |i(iV) 2 - iF>» [A„ A v \ - \{D»A»f 
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where S^, s is the five dimensional Chern-Simons action 
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S 5 CS [R] = ^2 / ^ Xap ^{R»d v R x d a R p + ^R^RxdaRp - *-R p R v R x R a R p ). 



(3) 



To derive equation (2), it is useful to employ 



ss cs [R] = g^2 / d5x ^ Xap Tr [3Rn (d v R x d a R p + i {R„R X: d a R p } - R v R x R a R p )} 



+ 



24tt 2 



d 4 xe^ Xa Tr 



<5i? M ( {R u ,d\R a } + -iR v R\R,j 



(4) 



For the special case where 5R P = d p 5u + % [R p , Sou] we have 



8u ( d p R v d x R a + -i {R p R v , d x R a } - -AR p d v R\R a 



(5) 



The resulting effective action takes the form 
Seff= &I Tr l D » UD * U ~ l ~ ^ - ^) 2 ] 

~24(U J^ Xap MJ,JMJ P ) 



+ 



J e fMs\a Trj-i?^ {R\, J a ] + R p R v R\J a 



48tt 2 J " [ 2 

+ -R^,JuR\J(j — JfiJuJ\R(j 



+ 



487T 2 

i 

48^2 



- 4i(L lja/ D' i U- 1 D v U + R^UD^U- 1 ) + 2{U~ 1 R tJLV UL^ - L pu R^)} 
+ 16[^ K,^] + i(ZV^) 2 + (A^) 2 )]} 



(6) 



where D^U = d^U + iR^U - iUL^ D^U' 1 = d^U' 1 + iL^U' 1 - iU^R^ J M = 
U^id^U and = -UJ^U' 1 . 

In order to make contact with the usual formulation of QCD, we write down 
the effective action for L il = R il = in eq. (6), to obtain 

Seff = MD.UD^U- 1 ) - — J e^Tr(J,JMJ p ) 

+ 4^2 / ^ Xa ^{^G^{G x ,J a } + (G^G iy + 1 -G^ + JpJ v )G x J a 
+ 4^2 / e^Tr\^(U- 1 G^ + G^){G x ,G^} + (G^G u -^G^ + G^)G x G u a 

+ 19P / Tr { 2 [MU-^V)] 2 - (D.U-'DM) 2 
-UG^ [D»U-\D V U] + 2{U- X G^UG^ - G^G^)} , 

(7) 

where G u = U^G^U - U'Hd^U = t/ -1 ^ + JjU. 



3. Stability analysis of classical solutions 

Our initial hope was to find nontrivial stable minima of the action (7). In 
the process of looking for such solutions, we observed some numerical instabilities, 
which caused us to suspect that the action (7) might be unbounded from below. 

In order to demonstrate that this is indeed the case, it is sufficient to show 
this for one trial set of functions, U tr i a i(x) and G tr i a i(x). In order to simplify the 
stability analysis, we will therefore begin with classical action without gauge fields, 

Gtrial(x) = 0. 

Had such an action lead to stable soliton solutions, it would amount to an 
(approximate) bosonization of free quarks in four dimensions, which would have 
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been an interesting result in its own right. As the following shows, this has not 
been attained in the present formalism. We shall comment later on what we believe 
might be the possible reasons for this. 

We therefore set the gauge field to zero, resulting in 

Seff [0,0, U\ = ^- 2 j Tr( J,J») - J e^Tr(J^J x J a J p ) 

i r (8) 

-192^ J / Tr[2( ^ J " )2 + (J ^ )2] 

Motivated by the Skyrme model, we are looking for a radially symmetric hedge- 
hog solution. We choose the classical solution to be a field of the form 

U c = e if{r)r - f 

where r are pauli matrices, the generators of some SU (2) subgroup of U (N c Nf) and 
/(r) is a radial shape function with boundary conditions /(0) = 7r, f(oo) = 0. The 
choice of the embedding will become relevant only if stable solutions exists, and 
then it should be discussed together with quantization of the collective coordinates. 
The Wess-Zumino term vanishes and the rest of the terms are given by 

(v-) 2 =(/" + ^-^) 2 

Next, we define an effective potential V e jj as minus the action divided by N c Nf 
and integrated over space only, 
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9A 2 r 

V eff (f) = —J dr [r\ff + 2 sin 2 /] 





oo 



+ 



48vr 



dr 



2 (rf" + 2f 



sin 



(2/)Y 



r / 



+ [r 2 (/ / ) 4 -4(/ / ) 2 sm 2 /] 



(10) 



Finding stable solutions of the action (8) is now reduced to functional mini- 
mization of the effective potential (10). 



Consider a family of trial functions 



/(r) = 2tan- i (-) 



(11) 



where a is a variational parameter which also determines the soliton size. The 
functions (11) satisfy the boundary conditions at both r = and r — > oo: 



f(r = 0) = vr 
f(r - oo) = 



;i2) 



The first condition is needed to ensure that the solution carries one unit of 
winding number, which in our normalization corresponds to one quark. 

For this family of trial functions, 



V eff (a) = 6A 2 a 



96a 



(13) 



Hence for a — > the potential is unbounded below. Thus, when attempting to solve 
the equations of motion, we will find that the soliton profile will be "squeezed" to 



zero width around the origin. This shrinking of the classical solution to zero size 
occurs despite the presence of both two- and four-derivative terms in the action 



A similar effect occurs already in the a model . Also there, the approximate 
action with up to four derivatives on the a and tt fields, constrained to a 2 + n 2 = f 2 
has its classical solutions "squeezed" to zero width, and with energy tending to 
— oo. 

In two dimensions, stabilization of the soliton is provided by the mass term. 
In four dimensions, one may add a mass term 



where vtlq is some scale related to the original quark mass in the QCD Lagrangian, 
(not quite the bare mass itself, as there is normal ordering to be performed; see 
Appendix of ref. [7]). Such a mass term does not provide the desired stabilization, 
however. For the trial function above, this term will have a divergent contribution 
to the potential coming from the integration over large r. Since the stabilization 
problem occurs at small distances, this large-r divergence due to the mass term 
cannot cure the problem. In order to isolate the large-r divergence, we will treat 
the problem of large r by putting a cut-off R. We expect that eventually such a 
cut-off will actually be provided by the confinement in QCD. 

Now the contribution of the mass term to the potential will be 



which tends to zero as "a" tends to zero, thus not changing the fact that the 



(8). 



2mS [1 - cos(/)] 



(14) 




(15) 
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potential is unbounded below for small scales a. 



4. Discussion and Interpretations 

There are various options to overcome these difficulties. The first is by choosing 
a different regularization scheme. This may change the coefficients of the dimen- 
sion four operators which appear in the Lagrangian. In particular it may change 
the terms in such a way that we will have a commutator squared as in the Skyrme 
model. This probably is the only known action which produces a positive Hamilto- 
nian. A second way out of this problem is to refer to non-perturbative corrections 
which will change the form of the coefficients in such a way as to get a Skyrme like 
action. 

We should remark, however, that in general we would not expect a scheme 
change to influence physical results, like the emergence of constituent quarks. It 
may happen, however, that due to the approximations made, we may be able to 
derive certain quantities in one scheme and not in another. 

Recall that in two dimensions, the scheme was completely fixed by requiring 
vector conservation and that the axial be the dual of the vector. The latter re- 
quirement was a result of our wish to have the bosonic version correspond to the 
fermionic one, and in the latter the axial is indeed the dual of the vector (see 
our work, ref. [8] for details). We do not have an analogous requirement in four 
dimensions as yet. 

Let us also remark that our classical configurations tend to be "squeezed" to 
zero size, and with energy tending to — oo. The troublesome part is at short 
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distances. But this is precisely the regime of high momenta, where our approxi- 
mations are inadequate, as we have neglected terms with six derivatives or more. 
So we either have to find a better approximation, or maybe exclude some short 
distance region. 

A final comment. We expect the effective action (7), after integrating out the 
gauge fields and taking trace over color, to yield an effective action in flavor space. 
But due to the non-positive nature of the potential that we discovered above, we 
do not expect, within the present approximation, to get the Skyrme model with 
the assumed standard positive-definite stabilizing four-derivative term 13 ' 41 . 
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